A ratchet driven by quasimonochromatic noise 
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The currents generated by noise-induced activation processes in a periodic potential are investigated 
analytically, by digital simulation and by performing analogue experiments. The noise is taken to 
be quasimonochromatic and the potential to be a smoothed sawtooth. Two analytic approaches are 
studied. The first involves a perturbative expansion in inverse powers of the frequency characterizing 
quasimonochromatic noise and the second is a direct numerical integration of the deterministic 
differential equations obtained in the limit of weak noise. These results, together with the digital 
and analogue experiments, show that the system does indeed give rise, in general, to a net transport 
of particles. All techniques also show that a current reversal exists for a particular value of the noise 
parameters. 



PACS numbers: 05.40. Ca, 05.10.Gg, 02.50.Ey 
I. INTRODUCTION 

The nature of microscopic engines, such as molecular 
motors, has been the subject of much research over the 
last five or six years. This recent activity was stimu- 
lated by the possibility of noise-induced currents [Q-Q, 
and was motivated to a large extent by the desire to 
model protein motors. These are proteins which are con- 
nected to a biopolymer and catalyze the conversion of 
adenosine triphosphate (ATP) to adenosine diphosphate 
(ADP). The energy released by this process is used by 
the motor protein to generate motion along the biopoly- 
mer in one particular direction. This is modelled as a 
microscopic object moving unidirectionally along a one- 
dimensional periodic structure It is this problem of 
rectifying processes at small scales that has stimulated 
most of the theoretical work in this area. A key ingredient 
is the presence of random Brownian forces. As a conse- 
quence it is natural to describe these stochastic ratchets, 
as they are frequently called, as a particle moving in a 
periodic potential subject to noise and to formulate this 
mathematically as a Langevin equation 



mx + ax + d x V(x, t) = £(t) , 



(1) 



where x(t) is the coordinate of the particle, a is a friction 
constant, V(x, t) is a periodic asymmetric potential and 
is the noise. 

Before discussing ([!]) in more detail, we should point 
out that there are at least three other reasons for the 
renewed interest in such systems. The first is a new gen- 
eration of experiments that can be performed in vitro || , 
which has led to the increased sophistication of the mod- 
els now studied. The second is the application of these 
ideas to non-biological situations at small scales — the 
realm of nanotechnology. Finally, this problem throws 
up fundamental questions concerning irreversibility and 
the second law of thermodynamics. These issues are dis- 
cussed extensively elsewhere H , so we will just make the 



essential point that it is only if detailed balance holds 
that we can use considerations based on the second law 
to decide that no coherent unidirectional motion is pos- 
sible §. 

The classification of different types of ratchets is in fact 
most easily carried out with reference to ([!]), since the 
terminology used to describe ratchets (correlation, flash- 
ing) is not always applied consistently. Nearly all studies 
neglect the inertial term in ([l]) and scale time by a, so 
that the coefficient of the x term is unity. Most studies 
have focussed on the cases (i) V(x,t) — V(x) is deter- 
ministic and is non- white (so that detailed balance 
does not hold), (ii) V(x, t) — V(x)<f)(t) where <j>[t) may be 
deterministic or random |Q and £(t) is white noise. Since 
the main prerequisite for any ratchet is that the system 
does not obey detailed balance and, since detailed bal- 
ance may be violated in many different ways, it is clear 
that many other forms are possible. 

In this paper we will study a ratchet of the type (i) 
discussed above, sometimes called a correlation ratchet. 
In section II the model is described in more detail: the 
potential is specified as is the type of noise — quasi- 



whose application to 
We study the model 



monochromatic noise (QMN) 
ratchets was first discussed in 
using analytic techniques in sections III and IV and by 
the digital and analogue simulations which are discussed 
in section V. We end with an analysis of results and con- 
clusions. Some of our preliminary work has already been 
reported M. 



II. THE MODEL 

In this section we will write down an explicit repre- 
sentation for a correlation ratchet acted upon by quasi- 
monochromatic noise. Having said this, it only remains 
to specify the potential V(x). We have already indicated 
that this function should be periodic and asymmetric. A 
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natural choice would therefore be a sawtooth potential: 
as shown in 0, a sawtooth potential is the one which 
maximises the current. However, in order to carry out 
analogue experiments we take the first few Fourier modes 
of such a potential leading to the form 

V(x) — 2 cos x + sin x + i cos(2.t) + i cos(3x) . (2) 

O -LU 

The sawtooth and the approximation to it, (||), are both 
shown in Figure 1, where we can see that we are mod- 
elling the sawtooth potential by rounding the corners. 
This not only makes it easier to reproduce in an analogue 
experiment, but it also removes possible singularities in 
the theoretical treatment: the "sharp corners" at the top 
and bottom of the potential would mean that the force 
that the particle would feel would be not well defined. 

Therefore the model is defined by the Langevin equa- 
tion 



x + V'{x) =£{t) 



(3) 



where V(x) is given by (||) and the noise is taken to 
be Gaussian with zero mean and correlation function 



(£H£(a/)) = 2D 2tt C(u>)S(u + uj') 



(4) 



We choose £ to be quasimonochromatic noise (QMN) 
since it exactly suits our purposes: the noise cannot be 
white, for reasons described in the last section and it 
has the physically appealing feature of having a peak 
at a non-zero frequency in its power spectrum (hence 
the name [|| QMN), while being simple enough to allow 
analytic progress to be made. Specifically, the noise is 
defined by 



c- 1 H = ( Ci ; 2 -^ 2 ) 2 + 4r 2 



(5) 



C{u>) is sharply peaked at the frequency (u>q — 2F 2 ) 1 / 2 w 
luq in the limit r « wj, and so we will frequently be 
working in this regime. This type of noise can also be 
viewed as the result of passing white noise through a 
harmonic oscillator filter: 



(6) 



(hence the name "harmonic" |10[|llfl , which is also often 
used) where the white noise r\ has strength D. 

We shall be mounting a three-pronged attack on the 
problem posed above: an analytic treatment based a 
small D approximation, direct digital simulation of the 
Langevin equation, to be discussed in section V, and an 
analogue simulation also to be discussed in section V. In 
the case of the analogue simulation, the first task, before 
simulating the Langevin equation itself, is to check the 
quality of the QMN produced by the analogue circuit. In 
order to accomplish this, we examine the QMN spectrum 
simulated by ffl). Two examples are shown in Figures 2 
and 3. 



It is evident that there is good agreement between the 
noise generated in the simulations and the theoretical re- 
sult given by (||). It is also clear from these figures that 
the shape of the spectrum changes considerably depend- 
ing on the values of the parameters T and u>. In fact, 
these figures illustrate the two regimes for QMN noise. 
The first (Figure 2) is an example of the case ujq > 2T 2 . 
The spectrum has a local minimum at uj = 0, rising to 
a maximum value at uj 2 = ujq — 2T 2 , and then falling off 
to zero as uj — > oo. As we have remarked already, in the 
limit ujq ^> 2r 2 the peak becomes narrower and better 
defined. Moreover, for values of uj such that w<wo, the 
spectrum is essentially flat and approximates well white 
noise. On the other hand Figure 3 illustrates the case 
uj 2 , < 2T 2 , where the spectrum has a local maximum at 
uj = and falls away to zero as w — * oo. So, in sum- 
mary, if the damping parameter T is small enough, the 
power spectrum has a peak at non-zero frequency. As 
r increases, the peak broadens and moves towards zero 
frequency. For T greater than a critical value of u>o/y2 
the maximum of the power spectrum is at zero frequency. 
Our aim is to see how the current changes as the noise 
parameters T and ujq vary. 



III. GENERAL FORMALISM 

In this section we discuss the approach we will use to 
explore analytically the generation of noise-induced cur- 
rents in the correlation ratchet introduced in the previous 
section. The method will involve an asymptotic analysis 
in the limit where the noise strength, D, tends to zero. 
To construct the asymptotic expansion it is first neces- 
sary to formulate the problem defined by the Langevin 
equation (^) in a different form. There are at least two 
different ways to proceed. One is to write down an equiv- 
alent Fokker-Planck equation. Since the noise £(i) in (||) 
is not white, it is first necessary to convert the process 
into an equivalent Markovian one with three degrees of 
freedom, (x, £, £), say. Thus the Fokker-Planck equation 
will have the form of time-dependent partial differential 
equation in three dimensions. We shall not pursue this 
method here, instead we will use the approach of ex- 
pressing the conditional probability (S(x — x(t)))ic as an 
average over all possible paths (or realizations of the pro- 
cess) with given initial conditions, denoted here by "IC" . 
These initial conditions specify not only the initial values 
of x(t), but also of £(i) and of £(t) at t — to. The explicit 
form for the path-integral is |12] 



P(x,t\lC,t ) = (5{x-x(t))) 



ic 



IC 



VxV[x]5{x - x(t)) 



(7) 
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where T>x is the appropriate measure defined so that P 
is correctly normalized and 



V[x] = J[x] exp-S[x)/D 



(8) 



Here S[x] is the action functional, which will be discussed 
in more detail below, and J[x] is the Jacobian of the 
transformation from r)(t) to x(t), for which we will not 
require an explicit form. 

The method for finding S[x] is discussed in some detail 
in [ jiljfl, but we can obtain it relatively quickly from (^) 
and (|6j) by first writing 



2r 

[x + V'(x)] + — [x + xV"{x)] + 
\ [x +xV"(x)+x 2 V"'{x)] = co 2 V {t) . 



(9) 



Since the noise 77 is Gaussian, white, with strength D, 
and has zero mean, the probability functional P[rj\ has 
the form exp — (1/4D) J dtrj 2 (t). A naive substitution of 
@ into this expression is sufficient to give the correct 
functional form for V[x] to leading order in D, namely 



where 

and 
A[x] 



V[x] = J[x] exp-.4[x]/A , 



A[x] = S[x]/<4 > A = D /^o 



(10) 



(11) 



dt 



2T 

+ V'(x)] + — [x + xV"{x)\ + 



[x +xV" \x) + x 2 V" \x)] 



(12) 



Having discussed the reformulation of the problem ex- 
pressed as the Langevin equation (||) as a Fokker-Planck 
equation or as a path-integral, we are now in a position to 
discuss the D — > asymptotics. In the case of the Fokker- 
Planck equation one may perform a WKB-like analysis, 
while in the case of the path-integral one may evaluate 
(0) by steepest descents, the paths which dominate the 
integral being those for which 



8A[x] 
Sx 







(13) 



From (|12|) it can be seen that a sixth-order nonlinear dif- 
ferential equation is obtained. The solutions of this equa- 
tion, subject to the appropriate boundary conditions are 
the instantons or optimal paths, x c (t), of the model. Sub- 
stituting this solution back into the action gives a number 
S = S[x c ] = 0JqA[x c ]. In either case, the WKB treatment 
or the steepest descent evaluation of the path-integral, an 
analysis of the conditional probability (ffl) leads to a rate 



of escape from one potential well to another which has 
the characteristic form A^exp— S/D, where AT will be 
termed the prefactor. In the ratchet we are interested 
in the current j which is proportional to the difference 
between the rates of escape from a particular potential 
well to the neighboring wells on the right and on the left. 
It is therefore reasonable that it should have the form 0] 



j = A [W+ exp(-S+/D) - N- cxp(-5_/L>)] 



(14) 



where the plus and minus symbols denote right and left- 
ward transitions respectively and A is the well spacing. 

In the next section we will calculate the actions in (14), 
by solving the sixth order equation obtained from (13]), 
numerically. However, in order to get some intuition for 
what may happen we will end this section by assuming 
that ujo is large (compared to the scale set by the cur- 
vature of the potential at the bottom of the wells) and 
obtaining the action as a power series in u>q 2 . In order 
to do this, we first rewrite the noise correlator (||) in the 
generic form 



C~V) = [1 + mr 2 uj 2 + k 2 tV 



where 



-2 1 



2r 2 
72 



K-2 



(15) 



(16) 



With the form ( |l5| ) , the action for a path, starting at the 
bottom of a well at x = a and ending at the top of an 
adjacent barrier at x — b, is given by |[14j] 



A = A[x c ] = / dx V + kit 2 / dx V'(V") 2 

J a J a 



,< 2 r j dxV [{V") 2 + v'v" 
fb 



- K ir 4 / dx {V'Y{V"'Y + 0(r b ) . (17) 

J a 

The 0(t 6 ) terms are also known, they are given in p4[ , 
and are proportional to k\ and k\K2- 

First of all, suppose that the spectrum is sharply 
peaked, uj^ 3> 2r 2 , then both k\ and K2 have magnitudes 
of order one. In this case the expansion ( |l7| ) is simply 
one in powers of r 2 = uJq 2 - The 0(t 2 ) corrections are 
necessary, otherwise A simply depends on the height of 
the barrier AV = J dxV'(x) and the asymmetry of the 
potential does not manifest itself. An interesting special 
case is when u>o is large, but Uq = 2r 2 . Then k\ is zero, 
and the second and fourth terms on the right hand side 
of ( p"7| ) vanish, as do all of the 0(t 6 ) terms. So in this 
case 



A = AV + t 4 dxV [{V") 2 + V'V" 



+ 0(r 8 ) [for 



2r 2 



(18) 
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If there were no 0(t 4 ) terms, it would be the case that 
when cuq 2 = 2T 2 , the action only depended on the barrier 
height, and so if the height of the barrier was the same 
to the right or to the left, we should not expect any net 
current. Moreover, if we plot the spectrum of the noise, 
for the particular value of ujq 2 = 2T 2 , we can see that 
it is very flat: the particle effectively is feeling a white 
noise which gives no ratchet effect at all. Since there are 
0(r 4 ) corrections in (18), this is not quite so however. 



In order to investigate this point in a little more detail, 
we have calculated the integrals in ( p"7j ) (and (|l^)) using 
V(x) given in (0). We find 



A, 



-± kit +a±K%T — a± k ± t + (J(t 



,(3) 



( 4 ) ^2 4 



where the coefficients a!± are given below: 



(19) 





a(D 


a< 2 ) 




a( 4 > 


+ 


4.62 


22.58 


681.41 


751.48 




4.62 


5.61 


59.06 


70.15 



Table 1. Numerical values of the integrals 
in @ for the potential (§). 

Let us focus on the particular values ojq = 10.5 and 
luq = 31.6 which we will use later. A short calculation 
using ([l9]) and the values of given in the table above 
shows, that to 0(t 6 ), A+ = A- when T = 6.84 (for 
ujq = 10.5) and T — 22.13 (for uiq = 31.6). If we had used 
the result ojq 2 = 2T 2 — valid for small r as indicated by 
( |l8| ) — the corresponding values of T would have been 
7.42 and 22.34 respectively. So we see that T = UJ0/V2 
is a reasonable estimate for the value at which A.+ and 
A- become equal when uj has the larger value (31.6), 
but it is considerably different in the case when ujq is 
smaller: uj — 10.5. We shall discuss the interpretation 
of the point where the actions for the left- and right- 
moving transitions become equal in more detail in the 
concluding section. 



IV. CALCULATION OF LEADING 
CONTRIBUTION 

In this section we calculate the leading small D con- 
tributions S± (or alternatively the leading small A con- 
tributions A±) to the current j in (^) for the case of the 
potential (|J). In the preceding section we illustrated the 
general idea by evaluating these actions for large u>o, but 



a general analytic treatment is not possible and we will 
have to resort to a numerical calculation of their values. 
^From (£L2f), the general form of the action functional is 



A[x] 



dt L(x, x, X] t) 



(20) 



The variation ( |T3| ) leads to an Euler-Lagrange equation 
of sixth order 



(-1) 



, 



(21) 



3=0 



where x^ = d 3 x/dP. A numerical solution will involve 
the decomposition of this equation into six first-order dif- 
ferential equations. A systematic procedure for achieving 
this is provided by the Hamiltonian formulation for the 
generalized mechanics given by (^o|) and (|2l| ) Jl5| . 

If we carry out this procedure starting from the La- 
grangian given by 



{[i + V'(x)} 

4 LJ. 



1 



2F 

— [x + xV"(x)] 



- 5 [x+xV"(x)+a?V m (x)]y 



(22) 



we find the following Hamiltonian 

H(x,p) = pix 2 + P2X3 + U0P3 2 - p 3 {u>l(x2 + V) + 

2r(x 3 + x 2 V") + x 3 V" + x 2 2 V'"} , (23) 

where x — (xi, x%, x 3 ) and p = (p\,p2,p 3 )- The action 
turns out to be fl5|| 



A 



pjdt 



Hamilton's equations have their usual form: 
. _ dH . OH 

Xi — -5 , Pi — ~~R 

dp l axi 



i = 1...3 



(24) 



(25) 



they are, by construction, equivalent to the Euler- 
Lagrange equations (|2l|). Using the Hamiltonian (|23|) , 
the six equations ( p5| ) yield 

±1 = x 2 (26) 

x 2 = x 3 (27) 
±3 = 2c4p 3 - {ul[x2 + V'{x x )\ + 2T[x 3 + x 2 V"{xi)\ 

+x 3 V"(x 1 )+X2 2 V"'(x 1 )} (28) 
pi - p 3 {uj 2 V"{ Xi ) + 2Tx 2 V'"(x 1 )) 

+x 3 V"'(x 1 )+X2 2 V""(x 1 )} (29) 

h = -Pi + p 3 {<4 + 2TV" (xi) + 2x2 V" (ai) } (30) 

p 3 = -P2+P 3 {2T + V"(x 1 )} . (31) 

For an escape problem, we are searching for solutions that 
provide the minimum of the action. Imposing the con- 
dition that the variation with time of the action is zero, 
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as before, H — 0, and following previous work GEL we 
choose as boundary conditions for the "uphill" solution 
(going from the bottom to the top of the potential) 

xi(-oo) = x min , ari(oo) = x max 
X2(±oo) = 

x 3 (±oo) = . (32) 

In order to solve (p6|)-(|3l|) in practice we have to truncate 
them to a large, but finite, time interval and use the 
boundary conditions 

Xi( Tj — x m i n ,xi(0) — 2T max 
x 2 {-T) = 0,x 2 (0) = 
x z {-T) =0,z 3 (0) =0 , (33) 

where we have used the time translation invariance of the 
equations. 

The TWPBVP subroutines developed for solving two- 
point boundary value problem were used |]l6| . In or- 
der to get convergence we used the following procedures. 
We linearized the equations (|2(])-([3l]) at the initial point, 
and changed the boundary conditions at that point, per- 
turbing them in the direction of the unstable manifold 
given by eigenvectors with eigenvalues having positive 
real part. We took as an initial guess a straight line join- 
ing the boundary points, and the solution of this modified 
problem was used as initial guess to the original problem, 
leading to a solution for the optimal path. 

Having found this solution, the action 

,o 

A= pjdt (34) 

J-T 

was calculated, and to minimize the effect of the cut-off 
effect we added the correction from integrating the linear 
expressions on the boundaries from — oo to —T and from 
to oo. When T is chosen to be large enough, the result 
is independent of its value. 

The results of this calculation will be discussed in sec- 
tion VI where they will be compared to the simulations 
that we have also carried out on this problem. We now 
turn to a discussion of these simulations. 

V. ANALOGUE AND NUMERICAL 
SIMULATIONS 

In this section we study our correlation ratchet using 
analogue and digital simulations. We have measured the 
variation of the current with T in two ways: analogue sim- 
ulation experiments on an electronic circuit and Monte- 
Carlo simulations on a digital computer. 

The first of these techniques |]rij-|l9|| involves the build- 
ing of an electronic circuit to model the system under 
study, the application of appropriate forces, and analysis 



of the response by means of a digital computer. The ab- 
sence of truncation errors makes analogue simulations es- 
pecially valuable for use e.g. with fast oscillating systems 
where the integration time (the time over which data are 
accumulated and perhaps ensemble-averaged) substan- 
tially exceeds the vibration period, as occurs with QMN. 
Digital techniques have the advantage that they can al- 
ways in principle be made more accurate than analogue 
methods, which typically achieve 2-3% accuracy, but the 
relative simplicity of analogue simulations and their real- 
ism (being much closer to a real experiment than a digital 
simulation) represent significant advantages. 

The electronic circuit used to model (||) and (||) is 
shown as a block diagram in Figure 4. The lower section 
is the harmonic oscillator used as a "filter" to convert 
H quasiwhite noise from a feedback shift-register noise- 
generator (20]]l|] into QMN. The QMN is then applied to 
the input of the upper part of the circuit, which models 
the ratchet potential itself. Although the basis of the cir- 
cuit is standard jl9| , several points of detail deserve am- 
plification. The force corresponding to the trigonometric 
potential of (^) is created using trigonometric identities 
to write it as: 

,m 29 . 4 . 6. 3 

V (x) = sm x + cos x sin x cos x H — sin x , 

10 3 5 

(35) 

so we can build the force using only two AD639 ICs [£lj . 
An inherent limitation of the AD639 IC is that it can 
only treat a restricted range of angles (±500°). To pre- 
vent its input from straying outside this range, provi- 
sion is made for resetting the circuit using DG303AC 
p2| switches (they have not been plotted in the block 
diagram) |17|-|l9[]. The voltage in the circuit correspond- 
ing to coordinate x was digitized with a 12-bit Microstar 
ADC fH, model DAP 3200a/415. Data analysis ex- 
ploited the on-board 100 MHz Intel 486DX co-processor, 
which was operated within a MatLab-based PC software 
system developed by Kaufman |24[] . 

The digital simulations were done using a specialized 
algorithm, described in ||, which we will briefly recall 
here. The particular structure of (||) and (||) poses, in 
principle, a problem if a simple minded algorithm is used 
in the integration: (^) is characterised by two time scales 
(uiq and r) and the integration time step (call it h) used 
in the digital simulations would be chosen in such a way 
that both u>Qh <C 1 and Th <C 1. Also, if r r is the typical 
relaxation time in (g), we should also satisfy h/r r -C 1. 
This latter inequality, for the typical parameters which 
are of physical interest, is normally satisfied as soon as 
the former one is: in other words, in a simple minded al- 
gorithm the constraint on the integration time step comes 
from (^J) rather than from (^), because the time scales in- 
volved by ujq and T are (much) smaller than r r . From 
the point of view of the problem we are trying to solve, 
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however, this would not be very efficient: we would be 
using most of the CPU time integrating the noise equa- 
tion (j^ rather than integrating the dynamical equation 
representing the model under study. 

The particular structure of (||) and (||) suggests that a 
specialized algorithm could be more efficient: the point 
is that (^|) is a linear filtering of an uncorrelated Gaus- 
sian noise. This means that the output of this equation 
(the variable £) is itself a Gaussian variable, of unknown 
intensity and correlation: hence, it may be possible to 
integrate (^) directly, working out the appropriate inte- 
gration algorithm, using the statistical properties of the 
Gaussian variable £. The algorithm used to integrate @ 
is the Heun algorithm [^5| which prescribes that we in- 
tegrate (H) with a couple of elementary steps, namely we 
first predict 



x(h) = x(0) + h[-V'(x{0))\ + r(h) 



(36) 



and then correct as 



x(h) = {x(h) + x(0) + h[-V'(x(h))] + r(h)}/2 . (37) 

In the expression above, we need to evaluate the quantity 
r{h) = Jq £(s) ds, which can be written || as 



where 



r{h) = A 31 £(0) + A 32 t(0) + w 3 



w 3 = B 31 zi + B 32 z 2 + B 33 z 3 



(38) 



(39) 



and where Zj are uncorrelated Gaussian deviates of av- 
erage zero and standard deviation one. Note that the 
quantity r(h) turns out to be a linear combination of 
Gaussian variables, as expected. Defining fl 2 = ujq — V 2 
and X± = — r ± \JT 2 — u>q, we have 



A 32 = 



2il 



1 e 



1 



A_ 



(40) 
(41) 



The expressions for B t j are very cumbersome, and we 
refer the reader to j|): note that J^] contains a misprint, 
the quantity 47rT/2f2 2 on the right hand side of Eq. (A14) 
should read TT/fl 2 . 

A warning is in order concerning the random noise gen- 
erator. The noise intensities of interest are fairly small 
compared to the barrier that the Brownian particle has 
to overcome to diffuse and generate a net current. It 
is then of great importance to make sure that the rare 
activation events are correctly generated, which implies 
that the noise generator should be particularly accurate 
in generating the tails of the distribution. The genera- 
tor used works by generating a Gaussian random vari- 
able, using the Ziggurath algorithm pq |, from flat ran- 
dom distributions obtained with a subtract and carry 
algorithm p7|,|28|]. 



The actions A± were calculated from the slope of plot- 
ting the logarithm of the mean escape time, calculated 
as the total time of observation divided by the number of 
transitions to the left or right, versus l/D. The current 
in the experiments is easily obtained by keeping track 
of the distance moved by the random walker and divid- 
ing it by the total simulation time. These results can be 
observed in the figures. 



VI. ANALYSIS OF RESULTS AND 
CONCLUSIONS 

In this section we wish to compare the theoretical pre- 
dictions of sections III, IV with the experiments. Our 
main aim is to understand the structure of the current 
( |l4| ) . This is made up of actions S± and prefactors A/± . 
For small D, the action dominates, so we begin by com- 
paring the actions calculated from (^) with the numeri- 
cal method discussed in section IV. For a typical value of 
LO{) = 10.5, the results are shown in Figure 5. This shows 
reasonable agreement between analytical and numerical 
results of solving the full set of equations of section IV. 
This fully justifies the approximation of section III, which 
is very useful given the difficulty of carrying out the nu- 
merical integrations in the method of section IV. 

Now we are in a position to compare the theoretical 
predictions with the digital simulation of section V. This 
is shown in Figure 6. The digital simulation shows the 
same trends as the analytical results, but there is a large 
amount of scatter. Nevertheless, the value of T at the 
point where A.+ = A- predicted by the digital simu- 
lations is in reasonable agreement with the theoretical 
value. 

The prefactors are, unfortunately, difficult to calculate. 
In fact, a calculation for QMN has not yet been carried 
out. The prefactors are, however, known for white noise 
and exponentially correlated noise (for small noise cor- 
relation time p9| ), and we use these in the expectation 
that they are a reasonable approximation to the true re- 
sult. In Figure 7, the current calculated from (|l4| ) using 
these prefactors together with an action calculated as 
in section IV, is plotted with the data from the digital 
simulation. A possible interpretation of the deviation of 
the theoretical from the experimental results could be the 
approximate prefactor. However it is clear that the expo- 
nentially correlated prefactor is an improvement over the 
white noise one, and this suggests that the correct QMN 
prefactor might give even better agreement. In any case, 
it is clearly demonstrated that a ratchet consisting of an 
asymmetric periodic potential plus quasimonochromatic 
noise forcing can indeed give rise to a net transport of 
particles. 

Table 2 displays the current obtained for loq — 10.5 in 
both the analogue experiment and the digital simulation. 



G 



In this case, although actions can be calculated analyti- 
cally, the parameter in the expansion for the prefactor is 
no longer small and one is unable to obtain an analytic 
expression for the current. The values luq = 10.5 and 
luq = 31.6 were chosen for technical reasons connected 
with carrying out the analogue experiment. 
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Table 2. Analogue (ja) and digital (jo) 
currents for ujq = 10.5 



There is no reason to assume that the current rever- 
sal (j = 0) necessarily occurs when A+ = A-, because 
the prefactors may cause some deviation from this lead- 
ing order result. However, from Figure 7 and Table 2 it 
seems that they do in fact occur at the same point — even 
though the magnitude of the analogue current is consis- 
tently less than that of the digital current. This may be 
because, as D — > 0, the action completely dominates or 
because the prefactors happen to be approximately equal 
at this point. 

In the literature, this problem has been already dis- 
cussed in M and |3(]]. In O it was found that there is 
indeed a net current in the system, and, working in the 
limit of small T/loq, the authors were able to show that 
the sign of the current changes as the curvature at lo = 
(i.e. C"(0)) changes with varying T. In J3(j the authors 
considered a model of the form (|J), but with an addi- 
tional white noise. However, as we have seen, there is 
no need to introduce an additional noise of this type in 
order to see a current reversal. The authors found that 
the current should change sign for at least two different 
values of T. 

Our aim has been to study how the current changed as 
the noise parameters T and ujq varied, for finite T/loq- We 
found that that the change in the spectral density curva- 
ture at lj = mentioned above is still the main effect in 
determining the current direction, in agreement with . 
We have included in the theoretical treatment higher or- 
der terms in T/luq: our result coincides with the result 
of in the appropriate limit, with a small shift in the 
transition point if T/luq is finite. The simulations which 



were carried out support the theoretical conclusions. We 
have not observed more than one current reversal experi- 
mentally (with fixed loq and varying T), but examination 
of shows that there is another solution (to 0(r 4 )) 
for which A+ = A-. It would be interesting to explore 
this regime in more detail experimentally. 
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FIG. 2. Noise power spectrum for u) — 31.6, V — 16.7. 
The jagged line is from experiment, and the smooth one from 
theory. The frequency / = uj/2-k. 



c 




2 4 6 8 10 12 14 

f (kHz) 



FIG. 3. Noise power spectrum for u>o — 31.6, V — 33.3. 
The jagged line is from experiment, and the smooth one from 
theory. The frequency / = cj/2-k. 



FIG. 1. (a) Sawtooth potential, (b) 'Smoothed sawtooth' 
potential. 
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FIG. 4. Experimental analogue circuit. 




FIG. 5. Analytic (small r) versus numerical (general r) 
actions for ujq — 10.5. Squares are for escapes to the left, tri- 
angles to the right. Symbols are from numerical integration, 
solid lines from analytic calculations. 
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FIG. 6. Digital simulations and theoretical actions, for 
uio = 31.6. Squares denotes escapes to the left, triangles to 
the right. Symbols are simulations, curves theory. 
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FIG. 7. Current for uj = 31.6 and A = D/u>o = 0.87. The 
squares are from digital simulation; the circles represent the- 
ory for a white noise prefactor; the triangles represent theory 
with the improved prefactor. 
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